function [ py ] = fuckfuckers( x, y, z, mesh, theta )
%FUCKFUCKERS Summary of this function goes here
%   Detailed explanation goes here

n = size(x, 1);
nTri = size(mesh.f, 1);

% F = zeros(3, 3, n);

py = zeros(3, 1);

for i = 1 : n
    
    p = [x(i), y(i), z(i)];
    
    for j = 1 : nTri        
        
        q0 = mesh.v(mesh.f(j, 1), :);
        q1 = mesh.v(mesh.f(j, 2), :);
        q2 = mesh.v(mesh.f(j, 3), :)
        
        T = triarea(q0, q1, q2);
              
%         qn = mesh.fn(j, :);
        
        v0 = fuckers(p, q0);
        v1 = fuckers(p, q1);
        v2 = fuckers(p, q2);
        
        temp = v2 * theta(3*j-2:3*j);
        
        temp';
        
        py(:) = py(:) + (v0 + v1 + v2) * T / 3.0 * theta(3*j-2:3*j);
        
    end
    
end

end

